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This document is an addendum to “One-point remapping of Lagrangian perturbation theory in the mildly 
non-linear regime of cosmic structure formation” (Leclercq et al., 2013). 


The remapping procedure, described in section II of 
the main paper (Leclercq et al., 2013), relies on the Eule- 
rian density contrast. As noted by previous authors (see 
in particular Neyrinck, 2013), in the Lagrangian repre¬ 
sentation of the large-scale structure, it is natural to use 
the divergence of the displacement field instead of the 
Eulerian density contrast S. This addendum provides 
additional comments on the one-point statistics of -0 and 
comparatively analyzes key features of ^ and 6. 

In the Lagrangian frame, the quantity of interest is 
not the position, but the displacement field T(q) which 
maps the initial comoving particle position q to its final 
comoving Eulerian position x (see e.g. Bouchet et al.^ 
1995 or Bernardeau et al., 2002 for overviews), 

x = q + 4'(q). (1) 

It is important to note that, though T(q) is a priori a full 
three-dimensional vector field, it is curl-free up to second 
order in Lagrangian perturbation theory (appendix D in 
Bernardeau, 1994 or Bernardeau et a/., 2002 for a re¬ 
view). We did not consider contributions beyond 2LPT. 
After publication of the main paper, Chan (2014) ana¬ 
lyzed the non-linear evolution of T, splitting it into its 
scalar and vector parts (the so-called “Helmholtz decom¬ 
position”). Looking at two-point statistics, he found that 
shell-crossing leads to a suppression of small-scale power 
in the scalar part, and, subdominantly, to the generation 
of a vector contribution. 

Let '0(q) = Vq • T(q) denote the divergence of the dis¬ 
placement field, where Vq is the divergence operator in 
Lagrangian coordinates, quantifies the angle-averaged 
spatial-stretching of the Lagrangian dark matter “sheet” 
in comoving coordinates (Neyrinck, 2013). Let 7 ^^,lpt 
and 7^^,Nbody be the one-point probability distribution 
functions for the divergence of the displacement field in 
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FIG. 1. Redshift-zero probability distribution function for the 
divergence of the displacement field computed from eight 
1024 Mpc//i-box simulations of 512^ particles. A quantitative 
analysis of the deviation from Gaussianity of these PDFs is 
given in table 1. The particle distribution is determined us¬ 
ing: a full N-body simulation (purple curve), the Zel’dovich 
approximation (ZA, light red curve) and second-order La¬ 
grangian perturbation theory (2LPT, light blue curve). The 
vertical line at ^ — 3 represents the collapse barrier about 
which values bob around after gravitational collapse. A 
bump at this value is visible with full gravity, but LPT is 
unable to reproduce this feature. This regime corresponds to 
virialized, overdense clusters. 


LPT and in Wbody fields, respectively. We denote by 
Vs the corresponding PDFs for the Eulerian density con¬ 
trast. 

In figure 1, we show the PDEs of for the ZA, 2LPT 
and full Wbody gravity. The most important feature of 
-0 is that, whatever the model for structure formation, 
the PDE exhibits reduced non-Gaussianity compared to 
the PDE for the density contrast S (see the upper panel 





2 


Model 

Vs 

P '0 

ZA 

2.36 ±0.01 

Skewness 71 

-0.0067 ±0.0001 

2LPT 

2.83 ±0.01 

-1.5750 ±0.0002 

iV-body 

5.14 ±0.05 

-0.4274 ±0.0001 

ZA 

Excess kurtosis 72 

9.95 ± 0.09 -2.2154 X 10“® ± 0.0003 

2LPT 

13.91 ±0.15 

3.544 ±0.0011 

N-body 62.60 ± 2.75 

-0.2778 ± 0.0004 


TABLE L Non-Gaussianity parameters (the skewness 71 and 
the excess kurtosis 72 ) of the redshift-zero probability dis¬ 
tribution functions Vs and V^ip of the density contrast S and 
the divergence of the displacement field '0, respectively. The 
confidence intervals given correspond to the l-cr standard de¬ 
viations among eight realizations. In all cases, 71 and 72 are 
reduced when measured from ^ instead of S. 

of figure 7 in Leclercq et al., 2013, for comparison). The 
main reason is that V 5 , unlike V^, is tied down to zero at 
S = — 1 . It is highly non-Gaussian in the final conditions, 
both in A/'-body simulations and in approximations to the 
true dynamics. For a quantitative analysis, we looked at 
the first and second-order non-Gaussianity statistics: the 
skewness 71 and the excess kurtosis 72 , 

7 i = ^ and 72 = ^ - 3, (2) 

where is the n-th moment about the mean and a 
is the standard deviation. We estimated 71 and 72 at 
redshift zero in our simulations, in the one-point statistics 
of the density contrast 6 and of the divergence of the 
displacement field 7 /^. The results are shown in table I. In 
all cases, we found that both 71 and 72 are much smaller 
when measured from V^p instead of V 5 . 

At linear order in Lagrangian perturbation theory (the 
ZeLdovich approximation), the divergence of the dis¬ 
placement field is proportional to the density contrast 
in the initial conditions, S{q), scaling with the negative 
growth factor, —Di{r): 

= Vq • = -Di(r)(5(q). (3) 

Since we take Gaussian initial conditions, the PDF for ^ 
is Gaussian at any time with the ZA. In full gravity, non¬ 
linear evolution slightly breaks Gaussianity. 7^^,Nbody is 
slightly skewed towards negative values while its mode 
gets shifted around ip ^ 1. Taking into account non-local 
effects, 2LPT tries to get closer to the shape observed in 
Wbody simulations, but the correct skewness is overshot 
and the PDF is exceedingly peaked. 

Figure 2 shows a slice of the divergence of the dis¬ 
placement field, measured at redshift zero for particles 
occupying a flat 512^-pixel Lagrangian sheet from one of 
our simulations. For comparison, see also the figures in 
Mohayaee et al. (2006); Pueblas & Scoccimarro (2009); 
Neyrinck (2013). We used the color scheme of the latter 


paper, suggesting a topographical analogy when working 
in Lagrangian coordinates. As structures take shape, ip 
departs from its initial value; it takes positive values in 
underdensities and negative values in overdensities. The 
shape of voids (the “mountains”) is found to be reason¬ 
ably similar in LPT and in the Wbody simulation. For 
this reason, the influence of late-time non-linear effects 
in voids is milder as compared to over dense structures, 
which makes them easier to relate to the initial condi¬ 
tions. However, in over dense regions where pj decreases, 
it is not allowed to take arbitrary values: where gravi¬ 
tational collapse occurs, “lakes” form and pj gets stuck 
around a collapse barrier, ip ^ —3. As expected, these 
“lakes”, corresponding to virialized clusters, can only be 
found in Wbody simulations, since LPT fails to accu¬ 
rately describe the highly non-linear physics involved. A 
small bump at t/; = — 3 is visible in 7^^,Nbody (see fig¬ 
ure 1). We checked that this bump gets more visible in 
higher mass-resolution simulations (200 Mpc//i box for 
256^ particles), where matter is more clustered. This 
means that part of the information about gravitational 
clustering can be found in the one-point statistics of pj. 
Of course, the complete description of halos requires to 
precisely account for the shape of the “lakes”, which can 
only be done via higher-order correlation functions. More 
generally, it is possible to use Lagrangian information in 
order to classify structures of the cosmic web. In par¬ 
ticular, DIVA (Lavaux & Wandelt, 2010) uses the shear 
of the displacement field and ORIGAMI (Falck, Neyrinck 
& Szalay, 2012 ) the number of phase-space folds. As 
pointed out by Falck & Neyrinck (2015), while these tech¬ 
niques cannot be straightforwardly used for the analysis 
of galaxy surveys, where we lack Lagrangian informa¬ 
tion, recently proposed techniques for physical inference 
of the initial conditions (Jasche & Wandelt, 2013; Jasche, 
Leclercq & Wandelt, 2015) should allow their use with 
observational data. 

Figure 3 shows two-dimensional histograms comparing 
Wbody simulations to the LPT realizations for the den¬ 
sity contrast 6 and the divergence of the displacement 
field pj. At this point, it is useful to note that a good 
mapping exists in the case where the relation shown is 
monotonic and the scatter is narrow. As pointed out by 
Sahni & Shandarin (1996) and Neyrinck (2013), matter 
in the substructure of 2LPT-voids has incorrect statisti¬ 
cal properties: there are overdense particles in the low 
density region of the 2LPT (5-scatter plot. This degen¬ 
eracy is also visible in the -0 > 0 region of the 2LPT 
V^-scatter plot. On average, the scatter is bigger with pj 
than with (5, in particular in overdensities (t/; < 0 ), since 
with LPT, particles do not cluster: p; takes any value 
between 2 and —3 where it should remain around —3. 

Summing up our discussions in this addendum, we an¬ 
alyzed the relative merits of the Lagrangian divergence 
of the displacement field and the Eulerian density con¬ 
trast 6 at the level of one-point statistics. The important 
differences are the following: 

1 . ^ being irrotational up to order two, its divergence 
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FIG. 2. Slices of the divergence of the displacement held, '0, on a Lagrangian sheet of 512^ particles from a 512^-particle 
simulation of box size 1024 Mpc/h, run to redshift zero. For clarity we show only a 200 Mpc/h region. Each pixel corresponds 
to a particle. The particle distribution is determined using respectively a full iV-body simulation, the Zel’dovich approximation 
(ZA) and second-order Lagrangian perturbation theory (2LPT). In the upper left panel, the density contrast S in the iV-body 
simulation is shown, after binning on a 512^-voxel grid. To guide the eye, some clusters and voids are identihed by yellow and 
purple dots, respectively. The “lakes”, Lagrangian regions that have collapsed to form halos, are only visible in the iV-body 
simulation, while the “mountains”, Lagrangian regions corresponding to cosmic voids, are well reproduced by LPT. 


contains nearly all information on the displace¬ 
ment held in one dimension, instead of three. The 
collapse barrier at t/; = — 3 is visible in 7^^ for N- 
body simulations but not for LPT. A part of the 
information about non-linear gravitational cluster¬ 
ing is therefore encoded in the one-point statistics 
of ^p. 

2. pj exhibits much fewer gravitationally-induced non- 
Gaussian features than 6 in the hnal conditions (hg- 
ure 1 and table I). 

3. However, the values of are more scattered than 
the values of S with respect to the true dynamics 
(hgure 3), meaning that an unambiguous mapping 
is more difficult. 

Note added. While this addendum was being refeered. 


the work of Neyrinck (2015) appeared. It uses the spher¬ 
ical collapse prescription for pj while checking various 
scales for the initial density held. The result is a fast 
scheme for producing approximate particle realizations. 
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FIG. 3. Left panel. Two-dimensional histograms comparing particle densities evolved with full A/’-body dynamics (the x-axis) 
to densities in the LPT-evolved particle distributions (the y-axis). The red lines show the ideal y = x locus. A turn-up 
at low densities is visible with 2LPT, meaning that some overdense regions are predicted where there should be deep voids. 
Right panel Same plot for the divergence of the displacement field '0. Negative corresponds to overdensities and positive '0 
correspond to underdensities. The dotted blue line shows the collapse barrier at -0 = — 3 where particle get clustered in full 
gravity. The scatter is bigger with '0 than with (5, in particular in overdensities, since with LPT, particles do not cluster. The 
turn-up at low densities with 2LPT, observed with the density contrast, is also visible with the divergence of the displacement 
field. 
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